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ABSTRACT 

We propose a novel approach to the problem of constraining cosmological initial con¬ 
ditions. Within the framework of effective held theory, we classify initial conditions 
in terms of boundary terms added to the effective action describing the cosmolog¬ 
ical evolution below Planckian energies. These boundary terms can be thought of 
as spacelike branes which may support extra instantaneous degrees of freedom and 
extra operators. Interactions and renormalization of these boundary terms allow 
us to apply to the boundary terms the held-theoretical requirement of naturalness, 
i.e. stability under radiative corrections. We apply this requirement to slow-roll 
inhation with non-adiabatic initial conditions, and to cyclic cosmology. This allows 
us to dehne in a precise sense when some of these models are hne-tuned. We also 
describe how to parametrize in a model-independent way non-Gaussian initial con¬ 
ditions; we show that in some cases they are both potentially observable and pass 
our naturalness requirement. 
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1 Introduction 


This old privilege notwithstanding, it took a long time to transform our gazing into the sky 
into physics. The study of the early universe truly left the mist of myth and speculation to 
become science only in the 1930’s, with the discovery of cosmic expansion. We had to wait 
until the 1960’s, with the discovery of the cosmic microwave background (CMB), to be able 
to discriminate between the Hot Big Bang and alternative cosmologies, and only in the early 
1990’s, with the detection of CMB’s inhomogeneities, did cosmology fully become a quantitative 
science. 

The celebrated WMAP survey P] has spectacularly conhrmed some general predictions of 
slow-roll inflation, and offered the possibility of signihcantly constraining alternative explana¬ 
tions for the primordial power spectrum. Future experiments may be able to go beyond the 
power spectrum and check for other features of the CMB, as, for instance, primordial non- 
Gaussianities. 

Present data already make it meaningful to ask about hner details of the mechanism that 
generates an almost scale-invariant power spectrum. We just mentioned one such detail: pos¬ 
sible non-Gaussian features. Another one is whether the correct initial state of the universe is 
the standard adiabatic “Bunch-Davies” vacuum [2j. Deviations from the standard inflationary 
vacuum offer the exciting possibility of observing very high-energy, “trans-Planckian” physics in 
the cosmic microwave background radiation, thanks to the enormous stretch in proper distance 
due to inflation [H]. This effect, which may present us with a real chance for probing string 
theory or any other model of quantum gravity, has received considerable attention, once the 
possibility was raised that these effects could be as large as H/M, with H the Hubble parameter 
during inflation, and M the scale of new physics (e.g. the string scale). 

Due to our ignorance of the ultimate theory governing high-energy physics, the most natural, 
model-independent approach to studying modihcations to the primordial power spectrum is 
effective held theory (EFT) jHE]. Using an EFT approach jl] concluded that the signature of 
any trans-Planckian modihcation of the standard inhationary power spectrum is 0(if^/M^), 
well beyond the reach of observation even in the most favorable scenario {H ~ 10^“^ GeV, M ~ 
10^6 GeV). 

What was absent from e.g. ref. |1] was a systematic EFT approach to initial conditions. 
Ref. |lj presented convincing arguments against the (in)famous a-vacua [Hj of de Sitter space, 
but it did not give a complete parametrization of hnite-energy, non-thermal states. That 
parametrization was given in mil, where the EFT approach was systematically extended to 
the choice of initial conditions. Refs. laiHi conclude that changes in the initial conditions for 
inhation are under control and may give 0{H/M) corrections to the primordial power spectrum. 
According to laiHi, these corrections are quite characteristic of UV modihcations to physics and 


can be distinguished from other corrections arising instead from IR changes in the vacuum state. 

In this paper, we shall argue that other constraints on the EFT of initial conditions make 
0{H/M) changes in the primordial power spectrum unnatural, in the same sense that a light 
Fermi scale is unnatural in the Standard Model. Before explaining further this point, we need 
to sketch the approach of refs. mu The most important difference between that and other 
approaches is that in Cl El initial conditions for modes of any wavelength are specihed at the same 
initial time t*. Other approaches give the initial conditions separately for each mode, at the time 
it crosses the horizon. The latter prescription is useful in the context of inflationary cosmology, 
but it obscures the held-theoretical meaning of the perturbation and/or initial condition: it does 
not easily account for the fact that after t* curvatures and energy densities are small, so the 
held theory is under control, and it does not easily translate into an EFT language. The former 
prescription, instead, leads naturally to a simple classihcation of initial conditions in terms of 
local operators dehned at the space-like boundary (i.e. initial surface) t = t*. It also allows one 
to rephrase the question of naturalness of initial conditions for inhation in terms of the usual 
held-theoretical notion of naturalness. 

In held theory we have a naturalness problem whenever the UV cutoh of the theory, M 
is much bigger than the observed value of the coefficients of unprotected relevant operators. 
For instance, a quadratically-divergent scalar mass term is unnatural (or, equivalently, 

hne tuned) whenever m -C M. In our EFT theory of initial conditions, we shall dehne a 
parallel notion. Initial conditions will be deemed unnatural whenever they require to hne-tune 
coefficients of relevant boundary operators to values much smaller than M^~^, where A is the 
dimension of the operator. 

Section 2 is devoted to summarize the EFT approach of mu- There, we also generalize 
their formalism to the case where 3-d boundary interactions are not added at an “initial” time, 
but they are used instead to parametrize an unknown period in the history of the universe. We 
shall do that to apply EFT methods to ekpyrotic/cyclic jH] or pre-big bang cosmologies HDl. 
In a sentence: we shall replace the unknown physics at the bounce -where the scale factor 
shrinks to zero to originate a spacelike singularity- with a spacelike (instantaneous) brane, 
whose world-volume supports local operators and additional degrees of freedom (see hg. [T]). 

Section 3 explains in greater details the concept of naturalness for initial conditions. Nat¬ 
uralness is then applied to constrain possible changes in the B.D. vacuum of inflation. They 
turn out to be either unnatural or IR universal; therefore, ill-suited to characterize signals of 
new high-energy physics. 

In section 4, we show that the spacelike brane that parametrizes the unknown Planckian 
physics at the bounce in cyclic/ekpyrotic/pre-big bang cosmologies can give raise to the correct 
power spectrum P ~ k~^, generically at the price of severe non-localities on the brane. We 



Figure 1: A schematic view of the EFT description of a cyclic cosmology. It plots the scale factor vs. 
time. The Planckian era around the bounce is replaced by an effective space-like brane (S-brane). 

show that typically these non-localities lead after the bounce to faster-than light propagation of 
signals originating in the pre-big bang phase. This pathology may not be lethal for these models, 
since Lorentz invariance is explicitly broken by the brane. Nevertheless, it is troublesome. It is 
avoided with a particular choice of interaction terms on the brane, which we proceed to show is 
unnatural according to our general prescription, since it requires the hne-tuning to zero of an 
unprotected relevant operator. The rest of section 4 is used to show how to write in local form a 
non-local boundary interaction on the brane, which can generate the correct post-bounce power 
spectrum. We show that the boundary term can be written in terms of local operators, but only 
at the price of introducing extra auxiliary helds that propagate only on the spacelike brane. We 
then conclude the section by showing that, besides leading to faster-than-light propagation, a 
generic boundary term also requires to hne-tune certain relevant operators, and it is thus also 
unnatural. Our argument conhrms and complements the analysis of m (see also pa)- 

Section 5 shows how to compute the effect of changes in the initial state of inflation that do 
not affect the power spectrum, but that do change the three-point function of scalar fluctuations. 
We show that the changes induced by cubic boundary operators can lead to observable non- 
Gaussianities, different from those studied in |Tni, without requiring undue hne-tunings. 

Section 6 studies the three point function of scalar fluctuations in cyclic cosmologies. Set- 






ting aside naturalness considerations, we show that, in this scenario, any pre-big bang non- 
Gaussianity is damped. Finally, we show that the particular form of the scalar potential used in 
cyclic cosmologies sets a bound on how close to the singularity is the space-like brane. Specif¬ 
ically, if we want to avoid uncontrollably large interactions, the pre-big bang cosmic evolution 
must be cut off and replaced by an effective brane at a time parametrically larger than M~^. 

In section 7 we summarize our hndings and conclude pointing out to possible developments 
of our formalism to other problems involving spacelike singularities. 


2 EFT of Cosmological Initial Conditions 

2.1 Inflation 


Differently from other approaches, here as in HI El we give initial conditions for modes of all 
wavelengths at the same initial time t*. 

The prescription starts by supplementing the EFT action describing all relevant low energy 
helds with a boundary term that encodes the standard thermal vacuum. To be concrete, we 
begin by working out the example of a massless scalar held in a time-dependent background. 
The 4d (bulk) action plus a 3d boundary term is 


^ = ^4 + ^ 3 , 


^4= ^ 


2 Jt* 


' r 


^.3 = 


d?x^h*{x) f d^y^Jh*{y)x*ix)K{x,y)xiy) 


( 1 ) 


Here h*^ is the induced metric on the surface t = t*. The role of S 3 is to specify the wave 


functional for the scalar y at t = G: 


^[x(a^)] = exp(i^3[y]). (2) 

Selecting an initial state for y means in this language to choose a particular k{x, y). For instance, 
in de Sitter space with line element 

1 


ds = a{ri) {—drj + dx^dx^), a{rj) = 


' Hy 


i = 1, 2, 3, —CX3 < p < 0, 


( 3 ) 


the standard thermal j2] vacuum is obtained by choosing 

I^IV 

Kik) = --—. 

1 -|- i\k\ri* 

Here, a tilde denotes the Fourier transform from space coordinates to co-moving momenta k 


( 4 ) 


= \/k ■ k) and rj* is the initial (conformal) time.^ This expression for k makes it clear that 
the choice of such initial time is conventional, since a change in rj* changes only k, not the wave 
functional. From now on, whenever needed, the standard vacuum functional will be called |0). 

^From now on, rj, rj* will denote the conformal time, t, t* will denote the synchronous proper time, and a() 
will always denote the scale factor. Notice that y* here is negative. 









2.2 Changing the Initial State 


Next, we want to find a convenient classification of changes to the initial state. This can be 
done by adding a new bonndary term to the action: S ^ S + AS 3 . To determine AS'3, we 
notice that, at any hnite time t after t*, we are insensitive to changes that only affect very low 
co-moving momenta k: co-moving momenta \k\ < H{t)a{t) correspond to pertnrbations with 
snperhorizon physical wavelength Ap > 1/H{t), which are nnobservable at time t. So, since we 
are interested in changes that can be observed in the CMB of the present epoch, we have an 
IR cntoff natnrally bnilt into the theory. This IR cntoff tells us that observable changes in the 
initial conditions can be parametrized by local operators: 


AR3 = E 

i 



( 5 ) 


Here O* are operators of scaling dimension Aj, M is the high-energy cutoff of the EFT and the 
/3ds are dimensionless parameters. The dimension Aj determines among other things how “blue” 
is the change in the power spectrum: the fractional change in the power spectrum is proportional 
to k^*~‘^. Since the EFT makes sense only for k < M, operators of high conformal dimension do 
not signihcantly change the observable spectrum. So, the most signihcant observable changes in 
the primordial fluctuation spectrum are parametrized by a few local operators of low conformal 
dimension. 

We just mentioned that the EFT needs a UV cutoff. This means that the operators Oi must 
be suitably regulated at short distance. In other words, they are local only up to the cutoff 
scale M. As a simple example, consider the dimension-four operator = {j3/M){dix)‘^- It has 
to be smeared at short distance, for instance by the replacement 


d.xdiX 9iX/(-9Va^(<*)Af^)9.X- 


( 6 ) 


Here f{x) is a smooth function obeying f{x) = 1, for x < 1 — e; /(x) = 0, for x > 1 + e; e is 
a small positive number. The scale factor a{t*) appears because we want to cutoff at M the 
physical momentum \k\/a{t*), not the co-moving momentum \k\. 


2.3 Power Spectrum 

As a hrst application, let us derive the change in the power spectrum of a minimally-coupled 
scalar held in de Sitter space, induced by the operator introduced in the previous subsec¬ 
tion [7j. The change in initial conditions AS'3, is equivalent to perturbing the Hamiltonian of 
the system by an instantaneous interaction Hj = —5{ri — r]*)AS^. So, the perturbed two-point 
correlation function is 



To first order in [3, the change is 

5G{k) = j d^xa{r]*){0\[O^{x), |x(fc, 0)|2]|0). (8) 

This quantity is easily computed in terms of commutators of free helds in de Sitter space, 
resulting in [ZllIlllIS] 

SG{k) = -^a{T]*)^lm [x+(p*, k)]‘^\k\‘^f{\k\ya\r]*)M'^). (9) 

The canonically normalized, positive frequency solution of the free-held equations of motion is 

X+(p, k) = - (1 + i\k\ri) exp{-i\k\ri), (10) 

\/w 

and (0||x(fc, 0)p|0) = H‘^/\k\^ is the unperturbed two-point function. For |A;||r 7 *| ~ 1, the effect 
of on the power spectrum, P{k) = {k/27r)^G{k), can be as large as 6 P/P ~ (3H/M, i.e. 
in the observable range when (3 is 0(1). For |A;||77*| > 1, the oscillating exponential gives a 
characteristic periodic signature jH]. The signal rapidly decays for |fc||?7*| 3> 1, so, it can be 
detected only if inhation lasts for a relatively short period; otherwise, O'^ would only affect 
unobservable super-horizon huctuation. 

2.4 Scalar Metric Fluctuations 

Of course, we are not really interested in the huctuations of a minimally coupled scalar. Those 
describe at best tensor perturbations of the metric. We want instead scalar perturbations of 
the metric. In the standard setting for slow-roll inhation, where the held content is the metric 
plus the inhaton <p. It is convenient to use the ADM formalism and decompose Qfj^y into 3-d 
metric shift and lapse N: 

ds^ = -N^df + hij{dx^ + N^dt){dP + N^dt). (11) 

Metric and scalar held split into background values plus huctuations as hij = exp(2p)(5jj + 7p), 
a{t) = exp(p), = 0(t) -f (p. The last equation means that the background scalar held, 0(t), is 
a function of the time t only. A particularly convenient gauge is 

hij = exp(2p -h 2 C){6ij -1- y^), dijij = 0, y** = 0, (p = 0. (12) 

In this gauge, the scalar huctuations are given by (, which is also the gauge invariant variable 
of ^n] (see also |I7j). The quadratic bulk action for ( is almost identical with that of a massless 
scalar held 

5 = 


IQtiG 

1 


dtd^x^ + eP{dCf 


p- 


(13) 


dgd x—a (p) + (^C) 


P 


IGttG 


(14) 











For future reference, we have written the action both in synchronous time t and in conformal 
time T]. Derivative w.r.t. t is denoted by a dot while a prime denotes derivation w.r.t. rj] G is 
the Newton constant. 

This formula shows that scalar fluctuations can be treated as a non-canonically normalized 
scalar helds. So, the equations derived in the previous subsections translate into formulas for 
changes in the CMB. One may worry that this conclusion has been reached with a very special 
choice of gauge. This is not the case, though. In the 4-d bulk, in any gauge, only the combination 
— {cf)/p)ip propagates. This combination is indeed the gauge-invariant dehnition of the Bardeen 
variable (see Hi and UTI section 10.3), and it is this quantity that appears in the action. On 
the initial-time boundary, instead, both C, and ip appear UTI, but the latter appears only as a 
non-dynamical held, without a genuine kinetic term. Schematically: 

Sboundary = J d^xF[({x,t*) , (15) 

where F[(,ip] is a quadratic function of its variables. Therefore, after imposing the non- 
dynamical (p equations of motion, dF/dp = 0, becomes a linear function of C, that can 
be plugged back into the boundary term to yield a quadratic function of C only. 

2.5 Cyclic Cosmologies 

We shall try to obtain an EFT description of cyclic cosmologies, independent of the details of 
their unknown high-curvature phase. In our description, the period of cosmic evolution when 
perturbations are generated, separates into three periods (see hg. 1 ): a slow contraction, where 
the scale factor evolves as airf) = {—ri/ri*Y, rj < —rj*, e -C 1^; a high-curvature bounce at 
\ri\ ^ r]*, the unknown dynamics whereof is replaced by a surface term (a spacelike brane, or 
S-brane); a FRW phase, typically a radiation-dominated phase with 0 ( 77 ) = 77 / 77 *, rj > 77 *. The 
scalar 0 is no longer a slow-roll inflaton. In the pre-bounce, pre-big bang phase, the potential 
is a very steep exponential over some range of 0 , for instance [TH] : 

V = —Voexp{—^JlQ7rG/p(j)), 0 < p <C 1, (16) 

and the background is 

j2p/8TrG 

a{r]) = p\ 0 (ry) = - log{-r]/r]*), (17) 

1 — p 

where the parameters Vq and 77 * are related by 77 * = \Jp{l — 3p)/87iGVo. After the bounce, 
the scalar decouples from the dynamics and one finds the action for scalar fluctuations in a 
radiation-dominated FRW cosmology. 


^In the pre-big bang model of cni, instead, e = 1/2. Here ry* is positive. 







To match the two phases, we need to add a boundary term, an S-brane, with nonzero tension. 
Its form is largely free; we can only say that it must contain at least the following terms: 


T+ 



d^xJh*_ + 


(18) 


Here T+ and T_ are constant tensions, needed to satisfy the junction conditions aX rj = rj* and 
7] = —r]*, respectively, is the determinant of the induced metric on the surface r] = ±r]* and 
... stands for gradient terms that vanish on the background C = 0- 

The nonzero tensions generate a stress-energy tensor, which is conserved but which violates 
the null energy condition (NEC). This is unavoidable if we want a spatially flat bounce. 

Thanks to the S-brane, we can have a bounce and we can write the action for scalar fluctu¬ 
ations in the gauge m as 


167rG^ 




r + {dcy 


+ 


P 



d^x 


/_^\ 2p/(i-p) 

V / 


J d^x{A[C{v*) -a{d)C{-V*)] + C{v*)F{d)Civ*)} ■ 


-C^ + {dCY 
(19) 


This is one of the most important equations in our paper, and it is worth of a few comments. 

First of all, the boundary term parametrizing the unknown high-energy physics at the bounce 
is made of two pieces. The hrst is the one containing the Lagrange multiplier A(a;). This is a 
non-dynamical held whose role is to link the value of the held ( after the bounce, at p = r^*, to 
its value before the bounce, at p = —rj*. ( can be continuous at the bounce, for a = 1, or it 
can jump, whenever a ^ 1. a{d) is a local function of spatial gradients. So, Civ*) is generically 
a function of ({—r]*) and its derivatives along the brane. Likewise, F{d) is a function of spatial 
gradients. Generically, it is nonlocal. Its role is to mimic super-horizon correlations induced 
by the unknown physics at the bounce. Superhorizon correlations do not necessarily signal 
acausalities, because our parametrization can also ht a “long” bounce, lasting for an arbitrarily 
long time. 

Notice that we did not introduce d(/dri in the boundary term. The reason is that d(/dri 
terms render the variation of eq. (uni) ill-dehned. So, whenever they do appear in a boundary 
term, they should be eliminated by an appropriate discontinuous held redehnition [Zj , or by using 
the bulk equations of motion to convert them into functions of ( and its spacelike gradients. 

We wrote eq. (Hi in a specihc gauge. To write it in a 3-d covariant form, we must re¬ 
express C in terms of the intrinsic curvature on surfaces dehned by = 0. To write it in a fully 
4-d covariant manner, we must furthermore express the position of the bounce in terms of a 
covariant (scalar) equation involving, say, 'ijj and the scalar curvature. 

We may worry about two features of this procedure. 






1 . A covariant equation for the position of the bounce reduces, in the gauge (IT^ . to 


r^ = G{C) + V*- ( 20 ) 

The function G is smooth and it vanishes at ^ = 0, but is otherwise unknown. So, 
generically, the bounce does not sit at a constant value of the time r] = rj*. This is 
not a problem since at quadratic order this bending of the brane only generates further 
quadratic terms of the form plus a hnite renormalization of the brane tensions T±.^ 
Explicitly, the induced metric on the brane is 


h*j = diGdjGgoo{G{x),x) + diGgoj{G{x),x) + djGgoi{G{x),x) + gij{G{x),x). ( 21 ) 
By substituting this formula into the universal term T f d^x\/hF we get, at linear order in 

c 


T J dW(3C + 3p'GcC), 






( 22 ) 


ri=G{x) 


This term induces a hnite renormalization of the brane tension, T —T(1 + p'G(^). At 
quadratic order in ^ we hnd as announced a term of the form C,PCi after eliminating time 
derivatives by the bulk equations of motion; 

T ( £xa^ [^2(9^0' + 2G^d^CN,a-^ + ?>{2p^ + p")GjC' + 6C' + 6GcCC'+ 


3p'GccC^ + ^(2C + 2p'GcC)^ 


(23) 


2. The functions a and F we introduced, are quite arbitrary. They can even break explicitly 
3-d rotations and translations, even though we will not do that in the following. So, 
neither the momentum constraint nor the Hamiltonian constraint impose any restriction 
on them. This is because our “initial time” brane is different from the boundary brane 
introduced in the context of de Sitter holography in ng. For us the brane is real. It 
either parametrizes a very specihc time in cosmic evolution (the bounce), or a specihc 
choice of quantum state at some time early in the inhationary epoch. Since we use the 
brane to specify an initial state, which is then used to compute expectation values of C 
at 7] —> cx) (that is today), we can change the value of p*, together with the form of the 
boundary terms, in such a way as to keep the late-time correlators constant. An example 
of this “renormalization group” is our eq. ©• So, our boundary terms at p* are the analog 
of arbitrary initial conditions for the RG group. In ref. ma instead, the late-time brane 
is a regulator, and the time evolution of all helds is hxed, because both late-time and 
early-time boundary conditions are hxed. The diherence between the two approaches is 
schematically shown in hg. 2. 


^The brane bending also introduces terms in G- As we have already mentioned, these terms must be canceled 
by held redehnitions, or expressed in terms of gradients via the equations of motion. 
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Figure 2: On the left, a schematic view of our brane, its position is represented by the vertical lines. 
By changing the operators on the brane and/or its position, we can change the initial conditions for 
the flow, hence the flow itself. On the left, the holographic brane of ref. m- Its boundary operators 
are constrained by requiring that the flow remains the same when its position is changed. 

3 Naturalness 

In particle physics, stability under radiative corrections is a powerful guide for constraining 
high-energy extensions of the standard model or of any EFT. By construction, an EFT has a 
built-in UV cutoff. In other words, the EFT only describes the low-energy sector of a theory 
whose UV completion is unknown. EFT is a powerful method whenever there is a large energy 
gap between a known low-energy sector (say the standard model) and an unknown high-energy 
sector -typically made of heavy particles- that decouples below a cutoff M. In this case, the 
effect of integrating out the high-energy sector is to introduce irrelevant operators in the EFT. 
They have dimension A > 4 and appear in the EFT with coefficients A change in 

the UV physics -as for instance a change in masses and couplings of the heavy sector- induces 
a small modification in the coefficients of these irrelevant operators. So, low-energy physics is 
shielded from changes in the UV physics, and EFT can be predictive even in the absence of a 
complete knowledge of high-energy physics. The problem lies in the relevant operators; those 










with dimension A < 4. Their coefficients must be much smaller than otherwise, the EFT 

would describe a trivial physics with no light states at all. A simple example of this pathology 
is a theory with two scalars, one heavy, with mass M, and another one with mass m. Below the 
energy scale M only one scalar propagates, and only if m -C M. Relevant operators are very 
sensitive to changes in the UV physics, the more so the lower their dimension. Generically, they 
appear in the EFT with coefficients Even if the tree-level values of these parameters 

are small, radiative corrections typically bring them to values Equivalently, any 

small change in masses, couplings etc. in the heavy sector bring these coefficients to their 
typical value. 

The best known naturalness problem in particle physics is that associated with the Fermi 
scale Mp 100 GeV or, equivalently, the Higgs mass. Suppose that the standard model holds 
up to a high energy scale M k, 10^® — 10^® GeV. Then, unless the Higgs mass is protected by 
symmetries (e.g. supersymmetry) it would be destabilized by radiative corrections, and driven 
to the UV cutoff. 

In our setting, we have a naturalness problem too. It arises because we change the cosmo¬ 
logical initial state by adding boundary terms to a bulk action as in eq. 0. These boundary 
terms are generically unstable under radiative corrections. So, even if we begin by modifying our 
theory by adding a (safe) irrelevant boundary operator^, we may end up generating dangerous 
relevant operators. 

3.1 Naturalness in Inflation I 

We begin by studying in detail the example given in section 2.3. There we computed the change 
induced in the fluctuation spectrum of a minimally-coupled massless scalar by the dimension-4 
operator 

0‘ = ^(a.xf. (24) 

This boundary term is not complete. It must be covariantized with respect to 3-d general 

coordinate transformations, which are not broken by a choice of initial conditions. The covari- 

antization is obvious and gives the following boundary term 

= ^ / d^x^/hh^^dixdjX- (25) 

Once we impose a shift symmetry y —> y + constant, is the lowest-dimension boundary 
operator quadratic in y. So, it is not unnatural to set to zero the coefficient of the relevant 
-and potentially dangerous- operator = aMx^- 

^Since the boundary is a 3-d field theory, all operators of dimension greater than 3 are irrelevant. 



In the gauge m, the boundary term (ESI) gives rise to interaction terms between y and (. 
In particular, we get a quartic interaction 

^ J d^Mv*)C‘^dixdiX- (26) 

This is a dangerous interaction. Its effect is best seen by writing the wave function at late time 
as a functional integral 

Xno^] = J (27) 

X,C are given free boundary conditions at the initial time rj*, while at late time, r] 0“, they 
are given hxed (Dirichlet) boundary conditions 

C(h)l»;=0 Cnowy X(h)|r?=0 Xnow- (28) 

This functional integral can be computed perturbatively either in terms of Feynman diagrams or 
using the Hamiltonian formalism of subsection 2.3. Results become clearer using the Feynman 
diagram approach. 

First of all, bulk interactions may be potentially dangerous, since naively the bulk action has 
interaction terms of the form {16nG)~^(d(d(, while the kinetic term is multiplied by the extra 
factor e = [see eq. (1141) ]. which is small in slow-roll inflation. Nevertheless, as shown 

in ESI, a field redefinition of ( eliminates all such terms and leaves only bulk terms of the form 

{IQirGyhXdCdC. 

The boundary interaction EHI), instead, does produce a dangerous effect. By integrating out 
X, this interaction generates a new boundary term for ( thanks to the diagram in £g. 3. 


\ 



Figure 3: Self-energy correction to the boundary “mass” terms. C’s are denoted by solid lines and 
x’s by broken lines. 
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The propagator for x is the standard one, computed in the BD vacuum 

ixi'ih -*:)) = ^7^3(1 + “ * 1 *^ 1 ’,') exp[!|A^|(),' - ),)|, ), > 


(30) 







By substituting this expression, computed ai r] = rj*, into eq. (123) we get 




A 

2M 


d^xa{T]*)C,'^ 


\k\ 


d^k 

<M (27r)3 2\k\^ 


\k\^{l + |fc| V') ^ I dV{v*)C". (31) 


Here we used a sharp momentum cutoff and the inequality M ^ H. Any other cutoff would 
give a like result. 

The effect of this term on the power spectrum is enormous, since the kinetic term for scalar 
fluctuations is multiplied by the small number e. If we canonically normalize ( by 


c = 



(32) 


we can use the same technique as in subsection 2.3 to get 


6P{k) 


Ca^(? 7 *)Im [v~^{7]*, /c)]^. 


k) 


H 


{l + i\k\vi) exp{-i\k\r]), 


\k\ 

a{r]*) 


<M. 


(33) 

With our normalization, C = PGM^/127re. So, for k ~ the power spectrum receives 

corrections 0{(5GM^/12'KHe). By imposing the observational constraint 6P{k)/P{k) < e we 
arrive at our hrst naturalness constraint on (3 


p < 127re\GM^ 


-A 

M' 


(34) 


Notice that the power corrections we computed in section 2.3 were at most 0{(3H/M). Since 
the term m is also linear in fd and as large as jdGM^/12'KHe, it dominates over the tree-level 
term whenever GM^/12'KeH > 1. Constraints of this kind were derived with a different method 
in da [HI. The method used here makes it clear that the bound comes from asking that the 
correction to the BD initial state is generic and robust under radiative corrections. Of course, as 
with any divergence, even that in eq. m can be canceled by appropriately choosing boundary 
counter-terms. The point is that this choice implies a hne-tuned UV completion of our EFT, 
which we should not assume without a valid reason such as a symmetry, or a better knowledge 
of the UV physics. 

Notice that the counter-term that cancels m is not the same that renormalizes the brane 
tension. The brane tension is renormalized because when we expand eq. (123 in powers of (C, we 
also get a linear term 

aA = ^ / d^xaivlCdiXdiX- 

This term can be canceled by changing the tension of the term 


(35) 












Equation m implies we can cancel the linear term dsni) by changing the tension as 


36T + = 0. 

487r2 

To cancel the quadratic term, instead, we would need to change the tension as 

96T + = 0. 

487r2 


(37) 


(38) 


3.2 Naturalness in Inflation II 


In the previous subsection, we introduced an extra scalar held, y, besides the inhaton. This 
was done to simplify our analysis and is by no means necessary. We could have introduced a 
boundary interaction involving only for instance® 


This boundary term can be covariantized by recalling that at linear order the 3-d Ricci curvature 
TZij is proportional to didj( plus a term linear in the transverse traceless huctuation 'jij. So, 
the 3-d covariant form of eq. (P|) becomes 


AR, = 


7 e 


d^x\/TiTZ. 


(40) 


M 8 irG 

Besides the quadratic term dnni), eq. (gni) produces, among others, a quartic interaction term 
In terms of the canonically normalized held v it reads 

dvrGy 


= 


Me 


[ Sx / v^d‘^v. 
J a{r]*) 


(41) 


The self-energy loop now produces a boundary term of the form 


971€ 


d^xa^{ri*)v‘^. 


Following the same steps we used to arrive to eq. (P|) we obtain a constraint on 7 


7 S 97re^(GM^) A. 


(42) 


(43) 


Covariantization of the boundary terms is just one of many ways in which dangerous bound¬ 
ary interactions appear. In the next subsection we will show that such interactions are induced 
at the one-loop level when we include the cubic vertices of the bulk theory. One interesting 

®The factor elSTrG)”^ is included for convenience so that for the canonically normalized field v the coefficient 
reduces to 7 /M. 










source of boundary interactions is the field redefinition needed to make all bulk self-interactions 
of ( O(e^). As it was pointed out in P2], all 0(e) bulk cubic interactions assemble into the form 



e 

SttG 


d7]d^xf{C) 




(44) 


where /(C) is a quadratic function of (. Since they vanish on the linearized equations of motion, 
in brackets, they can be canceled by a local held redehnition of the form C C + /(C) [ 12 ] • 
When substituted in eq. (Ei, this held redehnition generates quartic boundary interactions 

J (45) 

The function /(C) is given in jT2], and we shall use its explicit form in the next subsection, to 
estimate the size of induced boundary interactions. Here, it suffices to notice that it contains, 
among others, terms like {a/2a'Y{dC,Y, which contain no factors of e. So, the interaction (1451) 
can generate a large boundary mass term through a diagram as in hgure 4. 

To sum up, small corrections by seemingly benign irrelevant boundary operators generically 
induce by radiative corrections large, dangerous, relevant boundary operators. This problem 
does not signal an outright inconsistency of non-BD initial conditions for inhation, but it makes 
their description in terms of an EFT unnatural, that is very sensitive to its UV completion. 

An equivalent way of stating the problem is simply that the hrst correction to be expected 
in a generic EFT of initial conditions is 


aM 


d^xa^{r]*)v‘^, 


(46) 


where a is a small coefficient at most of order eH/M (because of experimental constraints!). This 
is the hrst universal correction to the BD vacuum that generically dominates over all others. Any 
UV modihcation of physics, be it strings, non-Lorentz invariant dispersion relations, or simply 
phase transitions in the late stages of inhation, reduces to this same term. In this perspective, 
modihcations of the primordial power spectrum seem ill suited to discriminate among diherent 
types of new high-energy physics. 

We could have guessed this result because (PF)|l is the lowest-dimension local operator that 
can be written with the held v and its derivatives. The only question is whether this operator 
is generated after all. Our explicit calculation answered in the affirmative. There is one last 
subtlety to explain. Locality depends on the variable we choose to parametrize scalar huctua- 
tions. The correct one is the canonical scalar held n, which has scaling dimension one. It is in 
terms of this held that one hnds that radiative corrections renormalize the coefficients of rele¬ 
vant operators as M‘^~^ etc. We can rewrite (14511 or any other local operator in v in manifestly 
covariant form using the 3-d metric hij. This expression need not be local in One possible 






covariantization of (gni) is 


aM 


SttG 




(47) 


where A is the covariant scalar Laplacian in 3-d. 


3.3 Naturalness in Inflation III 

The strongest naturalness constraints on the coefficients of the boundary irrelevant operator in 
(eq. IHH]) arise when we take into account the boundary nonlinear self-interactions of the held C,. 
In this section we show that the constraints obtained in this way are the same as those obtained 
in ref. OS from backreaction considerations. 

In [12], the cubic interaction terms of the held C, were found to be 

e 


= 


SttG 


d'^xd'qa (v) eCiC) + /(C)°C + 


Here is the term inside brackets in eq. 

/(C) = <^ + - 


iv) 


CC + 


,-2 


iv) 


(dcr, 


(48) 


(49) 


and we have omitted terms of higher order in e. These cubic interactions give UV divergent 
Feynman diagrams that induce new marginal and relevant boundary operators of the form 

^ / -.ri O / ^ ^ ^ ^ 


= J- 


dW(r/*)CC', 


= aM- 


d'^xa'^ {■)]*)( 


(50) 


SnG J SnG 

The operator must be interpreted as explained in section 2.5: the term on the boundary 
is incompatible with a consistent variational principle, so it has to be eliminated with a held 
redehnition. Equivalently , if we work only to linear order in 7 , it can be transformed into a 
term of the type by using the unperturbed BD boundary condition, Civ*) = ^C(h*)- 

The same procedure used in section 2.3 shows that these boundary operators lead to modi- 
hcations of the primordial power spectrum of the form 

SP‘P(k) M 


■yg^^^kr]*), and 


P{k) 


~ a —q 


( 1 ) 


{kr]*), 


(51) 


P{k) 

where and g^^'^ are both 0 ( 1 ).® 

The clearest way to exhibit the renormalization ehect of the cubic interaction on the bound¬ 
ary Lagrangian is to perform a quadratic held redehnition on (, of the form C C + /(C) !13j . 
This redehnition eliminates from the bulk Lagrangian the interaction proportional to /(C), leav¬ 
ing only cubic terms of O(e^) in the bulk. On the boundary, on the other hand, the ehect of 
the held redehnition is to produce new cubic and quartic interactions: 


d^xa{ri*)(d‘^( 


d^xa{g*) [C + /(C)] [C + /(C)] 


^Explicitly: = cosy — y ^ siny, y(^^(y) = y ^[2cosy — (1 — y^)y ^sin2y]. 


(52) 













Using the explicit form of /(C), eq. (liHll . we see that new boundary action contains, among 
others, the terms 


AS’l 


7 g 
MSnG 


d^x-Cd^ (CCO , 




7 g 
MSnG 


d'^x- 




M'd'HdCf 


(53) 


These interactions induce boundary operators of the form (EOD, through the Feynman diagram 
shown in £g. ^ 



Figure 4: Effective boundary two-point vertex induced at one loop by the field redefinition C ^ C+/(C)- 
The quartic interaction is proportional to 


By plugging the interaction AS^ in the vertex of this diagram we arrive at a contribution to 
the effective action for C: 
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MJ (27r)3 
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k 7] 


(54) 

(55) 

(56) 


where in the third line we have considered only the most divergent term, and we used for ( the 
scalar held propagator given in eq. PH). up to an normalization factor SnGe 

{({k, v)C{-k, v')) = ^^{x{k, v)x{-k, y'))- (57) 

e 

The induced boundary term is the operator in eq. dSOD, with a coefficient 


a 


SnGM^ 

rsj - 

^ H 


(58) 


A similar calculation shows that, by plugging the interaction in (jSSD in the vertex of 
the diagram in fig. El one generates the operator 0^°^ given in eq. (jSDI), with a coefficient 

7 SnGM^ 


7 




(59) 


By asking that 5P^^\k)/P{k) and 5P^^\k)/P{k) are within the presently acceptable deviation 
from scale invariance in the power spectrum, i.e. that they are smaller that the combination of 















slow-roll parameters (e + rj), we get, from eq. (15111 . the constraints fonnd in usi for the coefficient 


7 : 




(60) 


Again, we want to stress that these are natnralness bonnds, in that they can be avoided by 
tnning the coefficient of the bonndary operators in eq. (jSOD to a mnch smaller valne than they 
would generically have in the presence of the boundary perturbation, eq. (jsni). 


3.4 Coda: Is the BD Vacuum Natural? 

In the previous subsections we applied standard techniques borrowed from held theory to study 
the naturalness of initial conditions that, below the cutoff energy M, differ from those specihed 
by the Bunch-Davis vacuum. One might wonder what happens in the absence of these modihca- 
tions: in particular, one might worry that the methods we used make even the unperturbed BD 
vacuum hne tuned, thus making the naturalness bounds we found less interesting. This would 
happen if one could hnd Feynman diagrams producing large boundary renormalizations that do 
not arise from the perturbation 7 , i.e. they do not vanish at 7 = 0. This is not the case, if one 
uses as the unperturbed vacuum, i.e. the one corresponding to 7 = 0 in eq. dSHD, the vacuum of 
the interacting theory. Indeed, the interacting BD vacuum is dehned by the functional integral 

^bd[C] = J K--] exp(i^), (61) 

where S = J drjL is the action of the interacting theory, and the integration in rj runs from 7 * 
to 7 = —cx) + ie. The Eucildean continuation in rj selects the BD (or Hartle-Hawking [ 20 ]) wave 
function. In the language of ref. [7] we are thus choosing “transparent” boundary conditions. 
Now, the wave function at any later time 7 ' > rj* is also dehned by the functional integral (161|) . 
with the range of integration ranging from g' to —00 + ie. In this representation, there is no 
boundary at rj*, so the held redehnition will not produce any large boundary term. 

By splitting the integration in two regions, before and after the “initial” time 7 *, this ob¬ 
vious result is reinterpereted as a cancelation between boundary terms arising from the held 
redehnition of the bulk action drjL, and the boundary terms dehning the BD vacuum. So, 
the only nonlinear interactions left will be bulk terms suppressed by powers of the slow-roll 
parameter. This result shows that the slow-roll inhation background is stable against both 
bulk and boundary radiative corrections. In the presence of the boundary perturbation dSHD, 
however, this argument breaks down, and new boundary terms containing f{() are generated 
by the held redehnition. 

Finally, we should point out again another important feature that distinguishes the BD 
vacuum from all others. As we have seen, a modihcation of the BD vacuum by new physics at 


a physical energy scale M generically manifests itself in the EFT by the boundary term (inD. 
Generically, one would expect the coefficient a to be 0(1), yet observation constrains it to be 
smaller than 0{eH/M). Since this constraint comes from the size of SP{k)/P{k) at \k\ = 
i.e. at \k\phys = \k\/a{ri*) = H, an alternative route to make the correction (HH]) compatible 
with experiment is to make inflation last long enough to stretch \k\phys above today’s horizon 
scale. This requires a{rjnow)/dij]*) > H/Hnow In other words, the change in the initial state 
of inflation behaves exactly as any other pre-inflation inhomogeneity. Moreover, its effect is 
generically described almost entirely by a single operator, eq. ini), no matter what originated 
the change. So, the BD vacuum is natural also thanks to the very property that dehnes inflation: 
wait long enough, and all changes will be diluted away. 


4 Naturalness in Cyclic Cosmologies 

4.1 Scale-Invariant Spectrum 


Ref. m convincingly showed that a scale-invariant spectrum of fluctuations cannot be produced 
before the bounce in the contracting phase described by eq. (ED with p -C 1. So, it must be 
created at the bounce. Both effects can be seen very clearly in our EFT formalism. 

We parametrized the unknown Planckian physics at the bounce in a cyclic cosmology in 
eq. (ED- By requiring that eq. (ED is stationary under arbitrary variations in ( and A, including 
those that do not vanish at the bounce, we get the equations of motion for ( as well as the 
junction conditions at the bounce 


(h, C" + 2-C' + - 0, 


rj* 


for 1 ] > 1 ]*, b{ri) = 


C{v,x) 

Kv) 

C(? 7 *,fc) = a(/c)C(-p*,/c), 


(62) 


ll Zl 

p \t]* 


-^C + paik)A 

dr] j 


= 0 . 


p/(i-p) 


for r] < T]*, 


= 0 , 


~V 


(63) 

(64) 


Before the bounce, the metric can be well approximated by a Minkowski metric, for p <C 1, and 
the general solution to the ( equations of motion is approximately a plane wave 


((ri, k) = 


(65) 


After the bounce, the general exact solution of the equations of motion is 


C{v,k) 


JL 

p* 




\ * 
P > P . 


( 66 ) 







(67) 


The wave function of the fluctuation ({t], k) is given by the functional integral 

^[C(h,^)] = J[dCdA] exp{iS[C, A]), 

where S is given in (unD. The boundary condition at r/ —>■ —oo is ({rj, k) —>■ Ae because the 

pre-big bang initial state is the standard Minkowski vacuum. We want to compute the visible 
primordial spectrum, so we must compute tI/[C(h) ^)] after the bounce, at 77 —> -|-cx). So, the 

value of C at the late-time boundary is hxed: C(h)^) C(^)- the quadratic approximation 
used in (US), the wave function is Gaussian 

^[CW] = nexp 

k 

The power spectrum is a function of the equal time, two point correlator of (: 

{((k)((k')) = + k')P(k) (69) 

mm')} = /i<iC]c(^c(*:')i>ticn7/K]i'i'ic]r = (2;r)3i’(A.+fc')jj^^. ( 70 ) 

To compute ■j{k) we perform the functional integral Since it is Gaussian, it reduces to 

= lim exp 

rj^+oo 
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C{v,k)C{v,-k) 


(71) 


7(/c) 


C\k) 


( 68 ) 


The boundary at rj ^ —00 has been discarded using the - 1 -ze prescription to select the positive- 
frequency Minkowski vacuum, that is by running the contour of integration in rj over Euclidean 
time in the far past. This prescription selects the true vacuum even in the interacting theory^. 
It is this prescription that generates a real part in 7 , which is naively purely imaginary. 

A major simplihcation occurs because we are interested in observable modes, whose wave¬ 
length is much larger than the Hubble radius immediately after the bounce: 27r/|fc| ^ 1 / 77 ( 7 *) = 
7 *. In this case one can use the boundary condition at late time, ( = ({k), to approximate the 
post-bounce evolution of ( as 


mk)^ak) + E'^+o(\k\r,), 

7 

2 vr 

7 ^ > 7 > h • 

\k\ 

(72) 

Before the bounce, boundary conditions at 7 —>■ — cx) set 



C( 7 , k) = A [1 - i\k\{r] + 7 *)] + 0(|/c| V), 

2 vr 

-TTl ^V<-V ■ 

\k\ 

(73) 


The matching conditions at the bounce, eqs. ()68I64I1 . are then easily solved to give 


E 


i\k\/p\a\^ + F 

1/7*-F-7|7|/7|«|2^^ ’ 


(74) 


^See ^21 for another application of this method to cosmology. 














By substituting into eq. (ini) we find 


7(/c) 


1 \k\/p\a\^ — iF 

SttG 1 — rj*F — i\k\ri* / p\a\^' 


(75) 


whence 


( 2 ^)! 


P(k) = 47rG 


(1 — p*ReF)^ + (p*ImF + \k\p* /p\a\^Y 


(76) 


k^ |A;|/p|aP — ImF 

Continuity at the bounce sets a{k) = 1, F{k) = 0. In this case we get P{k) oc |/cp + ? 7 *|A:|^ ~ 
instead of the scale invariant spectrum P{k) oc cst. This result agrees with ref. jllj . 

This means that all the features of the primordial spectrum are due to the unknown physics 
at the bounce, where the NEC is violated. This is an unavoidable CH, unpleasant feature of 
cyclic cosmologies, that makes the whole pre-big bang phase almost useless. Almost, but not 
completely. We can consider p* as “the earliest” time and integrate out all the pre-big bang 
evolution, but in doing so, we would need nonlocal operators on the brane to generate a scale 
invariant spectrum. The long pre-big bang phase, instead, allows for one good choice of local 
coefficients on the brane: 


F{k) = 0, a{k) = 


2^*2 


a\k\‘^p 

Vp 


(77) 


where the constant a can be related to the inflationary parameters H and e introduced in section 
3 by comparing eq. (Hi to the standard inflationary power spectrum, P{k) = 47rGi7^/(27r)^e, 
which leads to ap* ^/e/H. 

Another natural choice is to require that ^ is continuous at the bounce. This sets a{k) = 
1, F{k) = {^/p<jp*‘^\k\)~^, with the same constant a as before. We shall call these junction 
conditions the “short bounce.” 

An interesting feature shared by all boundary conditions is that C, keeps evolving well into the 
FRW phase, because the coefficient E of the decaying part of C must be much bigger than (^{k) 
to create the correct scale-invariant power spectrum. For the junction conditions in eq. dZZI), 
for instance, 

(78) 


This evolution could be in itself a problem for the cyclic cosmology scenario. Eq. dZHl) shows 
that right after the bounce, in the radiation-dominated FRW phase, fluctuations are minuscule, 
compared to today. So, they can be affected by tiny inhomogeneities developing in the FRW 
phase, which we have ignored in our formalism. This is another potential source of fine-tuning, 
besides that due to our ignorance of the physics at the bounce, but we will not investigate it 
here. 

Finally, the post-bounce evolution effectively erases any non-Gaussianity that could exist at 
the bounce. We shall discuss this effect in more details in section 6. 






4.2 Causality Concerns 


The parametrization of the unknown Planckian physics at the bounce given in eq. (HD contains 
instantaneous interactions at the bounce time ± 77 *®, so acausal propagation of signals may 
occur. To check if this pathology does occur, consider a plane wave propagating before the 
bounce along the direction x. The metric before the bounce is slowly contracting when p -C 1 
[see eq. (ED] so we approximate it by a Minkowski metric as in the previous subsection. The 
wave is then 

= >)<->)•■ (79) 

Before the bounce causality, that is propagation inside the light cone, means ({r],x) = 0 for 
X > 7] + 7]* + constant. We can set the constant to zero with a translation of the coordinate x. 
Thus, causality implies that /(a;) is analytic in the upper half-plane Imo; > 0. 

After the bounce, the general exact solution of the equations of motion is (jnSD- 
Before the bounce, for the wave given in eq. dZHD, 00 = k. We choose Q;(a;) real at Imo; = 0. 
In this case, the junction conditions reduce to two linear equations for C, D 

(1 + 7\k\)c + (1 - 7\k\)D - F{k){C + D)- 7—^{C + D) = 0, C + D = a{k)f{k). (80) 

pa^yk) 

The coefficient of the right-moving wave is C for k > 0, and D for k < 0. Since k = 00 , 
this means that such coefficient always multiplies the term t] exp{—iuri). Call it E{ui). /,From 
eq. dHOD we hnd 

E{uj) = e{u)C{uj) + e{-uj)D{u) = fl-icj- E((j) - f(cj). (81) 

2u y pa^yuj) J 

Whenever E{uj), a{uj) are analytic in the upper half-plane Imo; > 0, and a{u) has no zeroes 
there, E{u}) is analytic in the same region, and no propagation outside the light cone occurs after 
the bounce. This happens for the matching conditions (HID. On the other hand, the “short” 
bounce condition, as well as most other acceptable ones, is non analytic: E{uj) = has 

a cut along the uj imaginary axis. 

The bounce explicitly breaks Lorentz invariance, so faster-than light propagation does not 
necessarily imply causality violations. After all, one needs a source before the big bang to get a 
signal that can propagate faster than light after the big bang. Yet, the pathology is troublesome. 
So, in the next subsection, we shall restrict out attention to the junction conditions ED and 
see whether they satisfy our naturalness criterion. Afterward, in section 4.4, we shall return to 
studying more general boundary conditions. 

^Identification of these two times implies an instantaneous bounce. Since we will find acausal signals propa¬ 
gating with arbitrarily high speed, this simplification will not affect our conclusions. 






4.3 Naturalness Concerns I 


The quadratic junction conditions given in subsection 4.1 need to be covariantized with respect 
to 3-d coordinate transformations. This is particularly simple for the conditions dZZD. It suffices 
to notice that at linear order d^C, ~ TZ. Then the action contains the term 
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d^xVhTZA 


(82) 


- 77 * 


By expanding this equation we get a coupling S'! ^ ~ The equal-time ( propagator is 

now oc so ^ produces a self-energy diagram as in fig. 3, which gives a one-loop correction 
of the form (recall that a{ri*) = 1 here ) 



(83) 


This correction destroys the scale invariance of the spectrum for all momenta |A:| < 0(VGM'^). 

Again, a scale invariant spectrum turns out to be possible by fine-tuning the operators 
inserted on the brane. Since they parametrize our ignorance of the high-energy physics occurring 
at the bounce, this means that a scale invariant spectrum depends on a very special completion 
of our EFT instead of being a robust, generic feature. 

Our analysis requires that we do not know in detail the history of the bounce. If we know 
what happens at the bounce, then we cannot invoke a genericity argument. To appreciate this 
point, imagine that we ignore the whole history of the universe before the bounce, i.e. we 
integrate out C,{ri) for r] < —rj* and we substitute for it the boundary term 


S3 
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IGttG 
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(84) 


A generic covariantization of this term would produce a self-energy correction resulting in a 
“mass” term ~ that changes the spectrum for all momenta \k\ < 0{GM^). In this case 

though, we do know the previous history of the universe. It is a long slowly contracting phase, in 
which the power spectrum l/\k\ is due to the masslessness of C, which is guaranteed by it being 
part of a gauge field: the metric. Essentially the same argument tells us that the scale-invariant 
power spectrum of inflation is natural. 


4.4 Naturalness Concerns II 


The function F{d) that we introduced in eq. (fn?]) parametrizes how the superhorizon modes get 
modified during the bounce. As we noted earlier in this section, in the case of a short bounce 
this function must be nonlocal to reproduce scale invariance at late times: 
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(85) 












The nonlocality of this interaction on the S-brane complicates the study of its stability under 
radiative corrections. 

Our strategy is to rewrite the part of the action containing F{d) as a fully local action, 
by introducing auxiliary helds localized on the S-brane. Covariantization of this local action 
introduces then certain universal interactions of ( with the additional helds, which allow us to 
study the radiative stability of this system in the same way as in the previous section. 

To introduce the technique, consider a three dimensional scalar with action: 

S[(p] = J (86) 

To render this action local, introduce an additional, massless held tf that couples to 0 as follows: 


5[0, (p] 


= / (VX 




(87) 


In terms of the variables 0, (p, this action is fully local. The held ip appears quadratically in it, 
so we can integrate p out exactly. After integration, we are left with an ehective action for 0 
only, which possesses a nonlocal propagator. A calculation of the diagram in hgure El gives: 




e//[ 


1 

47r 


d?x4>- 


A' 


+ ..., 


( 88 ) 


The ellipses denote terms of order 0^ or higher, that we do not need for computing the power 
spectrum. The quadratic term is indeed the nonlocal action of eq. dHSl)- We note that this 
mechanism depends crucially on the masslessness of the scalar p. If radiative corrections induce 
a mass m for (p, the form of the action eq. (I86|l is spoiled for all momenta \k\ < m. 



Figure 5: The one-loop, UV finite diagram giving rise to the effective action in equation (IS8I) . The 
solid line denotes p, whereas dashed lines denote 0. 


In our case, the fundamental scalar is the (^-held. We cannot exactly use the above mechanism 
with 0 —because the action in eq. (HZl) cannot be covariantized consistently. Indeed, the 
interaction term (p'^ is uniquely covariantized by the local expression y/hp'^. At zeroth order in 
the metric fluctuation however, this term gives a mass to p, which makes the one-loop induced 
term regular in the infrared, instead of being divergent as \k\~^. 

For this reason we must resort to a more complicated system to generate an action of the 
form (lHF)|l for (. The following covariant action satishes our requirements: 


s[C,p,i^] 


d^xy/h -|-00^A(y9)0-|-0A0 


(89) 










Here is a scalar, -0 is a spinor, and A is the covariant Laplacian. Notice that these helds have 
non-minimal kinetic terms. To see how the effective action for ( looks like we proceed in two 
steps, by first integrating out xjj and next (p. Integrating out ^|J at one-loop, as shown in figure IHl 
we arrive at an effective action for (p and C, that to quadratic order reads: 


Next, integrating out ip the effective quadratic action for ( reduces to: 

By appropriately choosing a and /3, we can make j 
at the S-brane is indeed (F{d)(. 


(90) 

(91) 

so that this interaction 


/ 

1 


\ 


\ 


/ 


Figure 6: The one-loop diagram giving rise to the effective action in equation (jnU). Dashed lines 
represent ip, whereas dashed-dotted lines represent p. 


The Lagrangian in dHi is our starting point for investigating if the induced action in (ED is 
stable under radiative corrections. The interactions in (P|) come from the covariantization of the 
action, with ( coupling universally to the other fields in the action. Now, we must check if any 
of these interactions feed down to dangerous relevant operators. For us a “dangerous operator” 
means an operator that spoils the scale invariance of the spectrum after the bounce, such as a 
large mass term for ip or ip. In principle, a mass smaller than the IR-cutoff, aH\now, is allowed, 
since we cannot probe the spectrum below this scale. After all, the function F{k) ~ should 
be regulated by that IR cutoff in order not to generate infinite backreaction at zero momentum. 
The point is that we will show that the natural induced mass is much larger than the IR cutoff. 

For example, we see that the covariant kinetic term for the '0-field in eq. (j89jl contains terms 
like: 


^ J d>xiia,c9‘i’- 


(92) 


This interaction gives rise to a mass term by the one-loop diagram shown in figure [71 Cutting 
off the momentum integral at the scale M, we obtain: 

'■ 




Ml 


Fxipip. 


(93) 


This relevant interaction dominates the original kinetic term for ip at momentum scales |A:| < 


Alp 











Figure 7: The one-loop diagram giving rise to the effective action in equation (1911) . Dashed lines 
represent "0 and the solid line represents (. 


More importantly, the one-loop induced operator in eq. PI|) is modified to: 

/\/f4 » 


(94) 


Now, by integrating out the field (p we obtain the following operator: 




1 

Mp 



(95) 


This term is clearly IR dominant over the original term in eq. (jnH), for all momenta \ k\ < M{M/MpY 
and it produces an unacceptably red spectrum after the bounce over the whole range of scales 
relevant for observation. To avoid the appearance of this term, we must tune the mass term 
of the fermion in eq. (0 to the much smaller value (by about thirty orders of magnitude!) 
discussed above, by adding a counterterm in the bare Lagrangian. This is the exact analogue 
of the mass instability of light fundamental scalars in the standard model. Our analysis here 
clearly indicates that the choice for the function F[d) in eq. is fine-tuned in an effective 
field theory perspective, and that only a very specific UV-completion of the short bounce can 
guarantee scale invariance at late times. 

One can show that similar considerations apply to other choices of boundary non-local oper¬ 
ators that reproduce the correct power-spectrum when inserted in eq. dzni). By an appropriate 
choice of auxiliary fields one can make these operators local. This always requires some of the 
auxiliary fields to have an unnaturally small value for the coefficient of an unprotected operator, 
that receives correction due to its universal interactions with the metric. 


5 Non-Gaussianities in Inflation 

5.1 Minimally Coupled Scalar in de Sitter Space 

A minimally coupled scalar models tensor fluctuations in the CMB, so it has an intrinsic interest 
besides offering a simple example of our method. The scalar action S '4 -|- S' 3 , given in eq. dH), is 
only the quadratic part of an action that can contain nonlinear terms both in the 4-d bulk and 
in the 3-d boundary. We shall consider here the effect of adding a boundary interaction 

ASf = j d:'xa?(r]')Xx^, 


(96) 





while keeping the bulk action S '4 quadratic. We could compute the effect of this term on late¬ 
time correlators in a Hamiltonian formalism, as in subsection 2.3, but we will use instead a 
Lagrangian approach, which is simpler for taking into account back-reaction effects and for 
application to cyclic cosmology. In this case, the effect at late time rj ^ 0~, is [cfr. eq. fITTIl ] 

nxnou.] = I ^ 97 ) 


The boundary condition at late time is x(0,a;) = Xnow{x) while at early time is obtained by 

3 

making the action S '4 -|- S '3 -|- stationary w.r.t. free variations of y, including those that do 

not vanish on the boundary. The resulting initial condition, written for the Fourier components 
X(0, fc), is 


x'iv*, k) + K{k)x{v*, k) + 3Aa(r/*) 


dH 


■xiv* J)x{v*, k - 1) = 0 . 


(98) 


In the bulk, y obeys a free equation of motion, whose general solution is 


x{ri, k) = A{1 — i\k\ri) exp{—i\k\ri) -|- 5(1 -|- i\k\ri) exp{i\k\ri). (99) 

3 

At tree level, the wave function & is \k[x] = exp( 2 S 4 -|- iSs + iAS^ ), computed on shell. 
We only need quadratic and cubic terms in the action, so we can expand the solution of the 
equations of motion as xiv^k) = Xoiv^k) -|- Xxi{'r],k), and keep only terms at most of linear 
order in A. The coefficients Xoiv^ k), Xiiv^ k) obey the boundary conditions 

Xo(0, A:) = x(/c), Xi( 0 ,A:) = 0 , x!o{.'n\k) + k{k)xo{,'n\k) = Q. ( 100 ) 

We need not write the boundary condition for Xi{'n*ik) because of the following reason. On 
shell and to linear order in A, we can integrate by part the free action and use the free bulk 
equations of motion to arrive at 

3 f 1 /* k" 

Si + Ss + ASi = \uia_^--J j^^a^{7])xQ{v,k)[xo{v,-k)+ ‘2\xi{v,-k)] + 

1 /* d^ k 

2 J k{k)xo{v*, k)][xoiv*, -k) + 2Xxiiv\ -k)] + 

J d^xa^{T]*)Xxl{v%x)Y ( 101 ) 

Boundary conditions (Unni) now make all terms within brackets vanish, except the hrst that 
reduces to x{~k)- Therefore, the action reduces to the standard quadratic term, which gives 
rise to Gaussian fluctuations, plus a cubic term obtained simply by writing the free held x ^-t 
7]* in terms of its late-time value x{k)- 


( 102 ) 







In this formula, we have dropped a real divergent term, which does not contribute to expectation 
values, since it disappears in taking the square norm of the wave function. Using eq. (Pl|l and 
the boundary conditions (Unni), the cubic term finally reads 

3 r d'^k ^ 

[xl{v*,k)] = X J j^a^{r]*){2nf6%ki + k2 + k3)Y[x{kj){l-i\kj\7]*)exp{i\kj\7]*). (103) 
The cubic non-Gaussianity is 


{x{h)x{k2)x{h)) = 

To first order in A this integral is 


I[dx]x{ki)x{h)x{h)\'^[x{kW 

I[dx{k)]\^[x{kW 


{x{ki)x{k2)x{k^)) = 6(27r)^5^(fci +/C 2 + fc3)A(fci, ^ 2 ,^ 3 ) n 




\{ki,k2,kz) = -2a'^(r7*)Im 


t=\ m 

3 

\{{l - i\kj\ri*) ex]){i\kj\ri* 
i=i 


(104) 


(105) 

(106) 


Notice that the oscillating exponent in A makes it vanish for momenta \k\ S> 1/|^7*|, i-e. for 
physical momenta \k\/a 3> H. So, the effect of the boundary term is significant only for 
momenta that crossed the horizon before the time p*. Physically, this effect can be understood 
by interpreting the boundary term AS^ as a “summary” of cosmic evolution prior to rj* -such 
as phase transitions, decoupling of heavy particles etc.. In this picture, any effect of this early 
history is washed away by inflation in all modes that keep evolving after rj*. This cutoff at rj* 
is another manifestation of the key property of inflating backgrounds, namely that they dilute 
away any primordial inhomogeneity. It also makes it difficult to detect such a non-Gaussianity 
unless Ha{r]*)/a{T]now) > Hnow 

For momenta \k\ 1/|?7*| and A real, eq. (j1 Objl becomes 

{x{ki)x{k2)x{k3)) = -{27if6^{ki + k2 + k^)^— ^ \ki\~^\km\~^. (107) 

^ l>m 


5.1.1 Estimating the Back-Reaction 

The coupling Ay^ has dimension 3 so it is marginal in the boundary EFT. Loops can only 
renormalize it logarithmically, so it can be naturally small. At O(A^), there exists a potential 
contribution to the boundary “mass” term / d^xof’x^. On dimensional grounds, we can easily 
estimate its coefficient as 0{X^M). So, as long as 

A^ < eR/M, (108) 


the induced “mass” term is not in conflict with experimental bounds. 








5.2 Scalar Fluctuations 


Scalar fluctuations ( have model independent non-Gaussianities usiini of the form 

(C^) = {27rfS^{ki + k2 + ^ ^ F{ki, ^ 2 , k^), (109) 

where F is a homogeneous function of degree —6 in the momenta. To compare with the results 
of the previous subsection, it is convenient to express the three point function in terms of the 
canonically normalized held v dehned in eq. (10211 . Then the size of (v^) is 0{V87rGeH'^k~^). 
Eventual non-Gaussianities in the initial conditions are instead 0{H^Xk~^). Glearly, for A large 
enough, the boundary non-Gaussianity can be dominant for momenta \k\ < 1/|?7*|. We must 
only check that the new cubic interactions do not introduce unacceptably large boundary “mass” 
renormalization. We already saw that this condition gives the bound A < ^JeH/M. Another 
bound comes from the mixed term due to the universal bulk cubic interaction UBI fivdvdv, 
where fi = 0{V87iGe). It is a dimension-5 operator that can induce at a boundary term 0{fi\) 
via the self-energy diagram shown in £g. 8 . On dimensional grounds the induced boundary 


/ \ 

I I 

\ / 



Figure 8: A mixed bulk-boundary interaction may generate a boundary mass term. The black square 
denotes the bulk interaction fivdvdv, the black circle denotes the boundary interaction Xv^. 

“mass” term is at most 0{XfiM‘^) J d^xa^x^- This gives another bound on A: 

A < M, ^ (8^G)-‘/l (110) 

For e ~ 10“^, M ~ 10^® GeV and H ~ 10^^ GeV this bound is even weaker than (ll() 8 |l . 

The ratio of the non-Gaussianity induced by boundary terms over the universal bulk term 
is XMp/y/eH . For the extremal value A ~ ^JeH/M, it becomes Mp/y/HM = O(IO^). So, there 
is dehnitely room for an observable signal here! A word of caution is necessary, tough. First 
of all, in order to see a clear signal we need that inhation does not last too long. Otherwise, 
thanks to the cutoff at \k\ ~ 1/|?7*|, the initial non-Gaussianity would only affect unobservable 
super-horizon huctuations. Moreover, in deriving the bound A < \J eH/M we only demanded 
that the induced change in the power spectrum is less than e. More stringent bounds on the 
power spectrum give a stronger bound on A ng. 











6 Non-Gaussianities in Cyclic Cosmologies 

6.1 Damping Effects 


As we mentioned earlier, in cyclic cosmologies, a scale invariant power spectrum for scalar 
fluctuations needs significant evolution after the bounce, well into the (radiation dominated) 
FRW phase. For the junction conditions d77|l . for instance, (now = ({k) = o'~‘^{\k\ri*)~^({ri*) = 
{\k\ri*)~^({—ri*). Suppose now that the wave function of the fluctuations just before the 
bounce contains a non-Gaussianity 


=l[exp 

k 


\k\ 

IQttG 


|C-(fc)r + Kki, k2, ksX-iksX-iksX-iks) 


C- = C(-h*)- (111) 


The three-point function can be computed as we did for inflation: 

,,3, ^ JWCt'l'IClP (16^ 

immc]? ~ w ' 


( 112 ) 


Assume for sake of example that bulk non-Gaussianities are negligible; then the late-time form 
of T is obtained by evolving ( with its quadratic equations of motion, and by expressing (^ in 
terms of (now- (- = (j\k\r]*(now 




A(fci, k2, k-i) 


Ilexp 


a^ri*^\kY 

IOttG 


|Ciioui(^)| T ^{kl) k2-! k^(now{k'i)(now{k^(now{k'i) 


aV I^i||^2||A:3|A(A: ^2, k^). 


(113) 

(114) 


Finally, the magnitude of {(now) is 




(115) 


So, any pre-bounce non-Gaussianity is damped by the extremely small factor {ari*\k\Y. Similar 
estimates hold for other choices of junction conditions. 

Evidently, only very dramatic non-Gaussianities have a chance to survive the smoothing 
effect of the bounce. In the next subsection we investigate whether this may happen. 


6.2 Non-Linear Interactions in the Slow-Contract ion Epoch 

There is one last subtlety about the size of non-Gaussianities in the cyclic models we are consid¬ 
ering that needs been addressed. It relates to the consistency of the perturbative expansion in 
this backgrounds. Looking at the scalar potential (ED, one may worry that for small p the non- 
linearities of the theory may become important, and the model may become strongly coupled. 
Indeed, if we consider the action to non-linear order, we may worry about terms containing 
self-interactions of the scalar, since on our background, for instance, VVp- 










Consider the situation in the harmonic (de Bonder) gauge, in which the kinetic term of 
the metric reduces to where h)) = h)) — and every component of the metric 

propagates. In this gauge, all interactions that involve the background metric are necessarily 
proportional to positive powers of H/Mp or (p'/Mp, since we are not solving for some of the helds 
in terms of other helds using the constraints. The only potentially large interaction is the cubic 
term in the expansion for the scalar held potential. Using the background solution, eq. (ED, 
this term reads 

-1 

The crucial fact is that this term is time-dependent. This means that we can trust perturbation 
theory before some critical time t]c that will be related to the value of p, but not afterwards®. In 
other words, we cannot follow the contracting solution arbitrarily close to the singularity. This 
means that, if we want to use perturbation theory, we are forced to place our S-brane at an 
earlier time than pc- 

An easy way to get an estimate for rjc is to consider the quartic coupling, This is a 

dimension four operator, and its dimensionless coupling is By requiring that this is 

smaller than 0(1), we arrive at the estimate pc ~ ■ Since we know that ri*a = 

and we need rj* > rjc our estimate of translates into p > {H/s/eMpYa"^ k, This should 

be combined with the requirement that the initial time is not outside of the effective held theory 
regime, p* > 1/M. 

7 Conclusions 

We presented a general method to describe unknown “dark” periods in cosmology, where high 
energy/Planckian curvature effects become important. These are, among others, the Planck- 
ian era at the beginning of inhation, and the singularity (bounce) in cyclic cosmologies. The 
unknown cosmic evolution during such times is replaced by a spacelike brane which can be 
endowed with extra degrees of freedom localized on it. These extra degrees of freedom allow 
for a general parametrization of the unknown UV physics by local operators of ever increasing 
dimension. 

There is an inhnite number of irrelevant operators, which parametrize the effect of the dark 
period and its unknown UV physics, as well as a hnite number of relevant operators. On the 
boundary, relevant operators have dimension strictly less than three. 

If the coefficients of the relevant operators are too small, they may be unstable against 
changes in the UV physics (or, equivalently, the details of the cosmic evolution during the 

^Recall that in the contracting solution we are considering, time runs from y = —oo to 77 = 0, at which the 
space-time is singular. 





dark period). The EFT parametrizes our ignorance of the dark period, so small coefficients 
of relevant operators signal a non-generic evolution, which is hne-tuned unless we have some 
additional knowledge of the dark period itself. From the point of view of the EFT then, such 
parameters are unnatural. 

We applied this naturalness requirement to constrain possible deviations of the initial state 
of inflation from the BD vacuum. We found that the generic signal of any such modification in 
the power spectrum of primordial fluctuations is a universal boundary “mass” term, no matter 
what is the origin of the UV modification. From this point of view, the power spectrum is not 
well suited to detect a clear, distinctive signal of new UV physics. For instance, trans-Planckian 
modifications to dispersion relations and late-time phase transitions during inflation would give 
an essentially identical signal. 

The naturalness requirement is quite powerful in constraining cyclic cosmologies. From this 
point of view, they are all hne-tuned. Many of them also allow for faster-than light propagation 
of signals and may be unstable against small inhomogeneities during the radiation-dominated 
phase of cosmic evolution. 

We also studied other signals of changes in initial conditions, besides those affecting the 
spectrum. We found that a short inhation still allows some initial non-Gaussianities to survive 
and be large enough for detection. Unlike the case of changes in the power spectrum, these 
non-Gaussianities can be large without generically inducing a large distortion of the power 
spectrum. 

Finally, we studied non-Gaussian features in cyclic cosmologies. We found that these cos¬ 
mologies generically damp all pre-bounce non-Gaussianities. We also showed that the request of 
no exceedingly large non-linear interactions in the pre-bounce, slow contraction epoch, demands 
that the S-brane cuts off cosmic evolution at a time parametrically larger than the Planck time. 

Our method extends the formalism set up in [2j to situations where cosmic evolution is 
known after and before a “dark” era, but not during it. The dark era is excised and substituted 
with an S-brane. It is fascinating to speculate that this method may shed new light on other 
cases where matter passes through (or ends on) a region that can be traded for an S-brane. 
The most interesting case is of course the singularity inside a non-extremal black hole, which is 
conjectured to be generically spacelike and of BKL [22j type. 

We want to end on a constructive note, by pointing out that naturalness could also become a 
guiding principle in the search for a UV completion of cosmic evolution near singularities, much 
in the same spirit as naturalness has guided particle physics in the search for a UV completion 
of the standard model. 
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